Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

2QQg  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2006  to  00-00-2006 

4.  TITLE  AND  SUBTITLE 

Sub-Bottom  Profiling  and  Geoacoustic  Inversion  Using  a  Ship-Towed 

Line  Array 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory ,4555  Overlook  Avenue 

SW, Washington, DC, 20375 

8.  PEREORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIEICATION  OE:  17.  LIMITATION  OE 

ARSTRAUT 

18.  NUMBER  19a.  NAME  OE 

OE  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

3 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Sub-Bottom  Profiling  and  Geoacoustic 
Inversion  Using  a  Ship-Toiued  Line  Array 

T.C.  Yang,  K.B.  Yoo,  and  L.T.  Fialkowski 
Acoustics  Division 

Introduction:  Bottom  properties  in  littoral  water 
can  be  very  complex  (i.e.,  not  horizontally  stratified) 
and  are  location  dependent,  as  indicated  by  core 
samples  and  seismic  surveys  taken  in  some  areas. 
Acoustic  signal  propagation  cannot  avoid  interacting 
with  the  bottom  in  shallow  water  areas.  Navy  opera¬ 
tions  in  littoral  regions  require  some  knowledge  of  the 
signal  propagation  loss,  which  determines  the  range 
that  a  (distant)  target  can  be  detected.  Knowing  the 
bottom,  thus,  becomes  critical  for  Navy  operations  in 
shallow  water.  Unfortunately,  there  exist  few  direct 
measurements  of  the  ocean  bottom  except  in  some 
isolated  locations,  since  core  sampling  and  seismic 
measurements  of  the  bottom  are  generally  expensive 
and  time  consuming.  Existing  data  bases  often  yield 
incorrect  predictions  of  propagation  loss  because  the 
data  entries  are  sparse  and  neglect  the  variation  of 
the  bottom  properties  along  the  propagation  path. 

A  recent  effort  at  NRL  is  devoted  to  developing  an 
acoustic  method  to  invert  for  the  bottom  properties 
(the  geoacoustic  inversion  method)  using  a  horizontal 
line  array  towed  behind  a  ship.  A  ship  can  cover  a 
wide  area  in  a  reasonable  time  (e.g.,  24  km  in  2  h).  The 
towed  array  inversion  method  can  provide  snapshots 
of  the  bottom  properties  under  the  ship  as  it  traverses 
an  area,  and  in  this  way,  build  a  data  base  for  the  Navy. 
This  article  summarizes  recent  advances  made  at  NRL 
and  experimental  verifications  of  the  method.^ 

Sub-Bottom  Profiling  and  Geoacoustic  Inver¬ 
sion:  Geoacoustic  inversion  correlates  the  signal  calcu¬ 
lated  using  an  assumed  bottom  model  with  the  acous¬ 
tic  data  received  on  an  array,  and  iteratively  alters  the 
parameters  of  the  bottom  (and  other  relevant  param¬ 
eters)  until  there  is  a  good  match  between  model  and 
data.  The  method  has  been  tested  widely  using  vertical 
line-array  data  and  works  when  there  is  some  prior 
knowledge  about  the  bottom,  such  as  the  number  of 
sediment  layers  in  the  bottom.  It  is  generally  known 
that  geoacoustic  inversions  can  have  ambiguous  solu¬ 
tions,  namely,  that  often  more  than  one  solution  can 
fit  the  data.  For  building  a  database,  the  inversion 
method  becomes  untrustworthy  if  multiple  data 
measurements  of  the  same  area  yield  different  bottom 
profiles.  Additional  constraints,  such  as  prior  informa¬ 
tion,  are  needed  to  narrow  down  the  inversion  results.^ 

In  the  seismic  community,  long  horizontal  arrays 
are  used  to  survey  the  bottom  using  acoustic  signals 
received  from  a  distant  source.  Focusing  on  bottom 


penetrating  acoustic  rays,  the  reduced  ray  travel  time 
is  deduced  as  a  function  of  the  ray  arrival  angle.  This 
relationship  can  be  used  to  determine  the  bottom  layer 
thickness  and  bottom  sound  speed.  In  the  acoustic 
community,  geoacoustic  inversion  is  preferred  to 
seismic  methods  for  two  reasons:  The  towed  sonar 
array  is  short  compared  with  the  seismic  array,  and  the 
source-receiver  range  is  fixed  at  a  few  hundred  meters. 
The  seismic  community  is  interested  in  very  deep 
bottom  structure,  whereas  a  shallow  (sub-)  bottom 
structure  is  more  critical  for  sonar  applications.  The 
geoacoustic  inversion  method  determines  not  only  the 
bottom  thickness  and  sound  speed  but  also  the  bottom 
density  and  attenuation  coefficients. 

NRL  developed  a  unified  signal  processing 
approach  that  combines  a  sub-bottom  profiling 
method  similar  to  that  used  by  the  seismic  community 
with  the  geoacoustic  inversion  method  used  by  the 
sonar  community.  By  beamforming  the  towed  array 
data,  the  beam  outputs  exhibit  an  angle  and  time 
(or  depth)  relationship  (Fig.  3),  from  which  one  can 
extract  a  reduce-time-angle  (or  slowness)  relationship 
similar  to  that  used  by  the  seismic  community.  Figure 
3  shows  a  direct  arrival  from  the  source  at  an  angle 
of  approximately  5°,  the  surface-reflected  arrival  at 
approximately  18°,  and  bottom  interacting  arrivals  at 
approximately  30°.  The  bottom  arrivals  show  a  return 
at  the  water-bottom  interface  and  a  return  from  the 
sub-bottom.  The  time  difference  is  related  to  the  sub¬ 
bottom  thickness.  Figure  4  displays  the  beam  intensity 
time  series  at  30°  for  data  collected  over  a  24-km  ship 
track;  this  shows  a  sub-bottom  layer  structure  (right 
figure)  in  good  agreement  with  that  obtained  from 
a  bottom  profiler  (left  figure).  The  arrival  angles  and 
bottom-layer  structure  provide  the  needed  inversion 
constraints  to  define  sufficiently  narrow  parameter 
search  windows,  yielding  an  inverted  bottom  profile 
for  water  depth  and  sediment  thickness,  bottom- 
sound  speed,  and  density  and  attenuation  (Fig.  5). 

The  inverted  result  at  the  edge  of  the  track  (designated 
as  site  4)  is  in  good  agreement  with  an  independent 
seismic  reflection  measurement  conducted  at  site  4. 

Summary:  A  method  of  geoacoustic  inversion 
incorporating  the  information  obtained  from  sub¬ 
bottom  profiling  is  developed  for  a  ship-towed  line 
array.  The  sub-bottom  profile  data  are  deduced  from 
the  acoustic  data  and  used  to  define  parameter  con¬ 
straints  used  for  geoacoustic  inversion  of  the  same 
data.  The  inversion  result  can  be  shown  to  be  ''unique" 
given  the  certainty  of  the  bottom  layer  structure.^ 

The  method  is  being  extended  to  use  the  ship's  own 
noise  as  the  acoustic  source.  In  this  case,  the  acoustic 
ray  arrival  angles  are  used  to  constrain  the  parameter 
search  windows  and  the  inversion  results. 
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FIGURE  3 

Beam-intensity  time  series 
at  different  beam  angles 
measured  from  the  forward 
endfire.  Effective  depth  is 
travel  time  x  1582  m/s. 
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FIGURE  4 

Effective  depth  determined  from  the  bottom  profiler  (left  figure).  Effective  depth  estimated  from  the  towed-array  data  superimposed  with 
the  profiler  data. 
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FIGURE  5 

Geoacoustic  inversion  results  for  the  water  and  sediment  depth  (upper  left),  sediment  and  bottom  sound  speed  (upper  right),  sediment 
and  bottom  density  (lower  left)  and  sediment  and  bottom  attenuation  (lower  right)  along  the  track  in  longitude.  The  numbers  on  the 
right  column  are  seismic  reflection  measurements  at  site  4. 
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